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Abstract 



It is proved in the paper that the non-conservative dissipative force and asymmetry of time reversal can be 
naturally introduced into classical statistical mechanics after retarded electromagnetic interaction between 
charged micro-particles is considered. In this way, the rational dynamic foundation for classical statistical 
physics can established and the revised Liouville's equation is obtained. The micro-canonical ensemble, the 
canonical ensemble, the distribution of near-independent subsystem, the distribution law of the Maxwell- 
5jT)' Boltzmann and the Maxwell distribution of velocities are achieved directly from the Liouville's equation 

without using the hypothesis of equal probability. The micro-canonical ensemble is considered unsuitable as 
the foundation of equivalent state theory again, for most of equivalent states of isolated systems are not the 
states with equal probabilities actually. The reversibility paradoxes in the processes of non-equivalent evolu- 
tions of macro-systems can be eliminated completely, and the united description of statistical mechanics of 
equivalent and non-equivalent states is reached. The revised BBKGY series equations and hydromechanics 
equations are reduced, the non-equivalent entropy of general systems is defined and the principle of entropy 
increment of the non-equivalent entropy is proved at last. 
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1. The fundamental problems existing in classical statistical physics 

Though classical statistical mechanics has been highly developed, its foundation has not yet been built 
up well up to now W. The first problem is about the rationality of the equal probability hypothesis or 
the micro-canonical ensemble hypothesiswhich is used as the foundation of equilibrium theory now. The 
hypothesis had got much criticism since it was put forward. In order to provide the hypothesis a rational 
• i— i ' base, Boltzmann raised the ergodic theory, proving that as long as a system was ergodic, the hypothesis 

of equal probability would be tenable. However, the study shows that the evolutions of systems can't be 
ergodic generally ^ . So the hypothesis of equal probability can only be regarded as a useful work principle 
Q_c without strict proof. As for the non-equilibrium statistical systems, we have no united and perfect theory 

at present. We do not know how a system to transform from a non-equilibrium state to an equilibrium 
one at present and how to define the non-equivalent entropy of general systems. Besides, there still exists a 
so-called reversibility paradox in the theory. Though lots of researches have been done, the really rational 
solution still remains to be explored. 

The reversibility paradox has been a long-standing problem. There exist two forms of the reversibility 
paradox at present. The first is the so-called Poincanre's recurrence. It was put forward by Zermelo in 1896 
based on a theorem provided by Poincanre in 1890 ' 3 - ) . According to this form, a conservative system in a 
limited space would return to the infinitely nearing neighbor reign of its initial state. The basic ideal of the 
proof is described as follows. For a conservative system, the Liouville's theorem is tenable, so the volume 
of phase space is unchanged in the evolution process. Because the volume of phase space is limited, the 
system should be recurrent after a long enough time's evolution. That is to say, the process is reversible. 
The second form was put forward by Loschmidt in 1876. Loschmidt thought that any micro-dynamic motion 
equation did not change under time reversal, or the motion of any single micro-particle was reversible, so 
after the velocities of all particles in a macro-system are reversed at the same time, the system would evolve 
along the completely opposite direction. Therefore, the process would be reversible. However, in the evolu- 
tion processes of isolated macro-systems, what can be observed is that the processes are always irreversible. 
Therefore, there exists the so-called reversibility paradox. Though a large numbers of explanations have 
been given up to now, none of them are satisfying. 



2. The Lorents retarded force and the basic hypothesis of classical statistical physics 



It is well-known that there exist two kinds of motion equations, one is for single particles, the other is 
for a statistical systems composed of a larger numbers of particles. The so-called micro-equation of motion 
that Loschmidt mentioned about in his time was actually the Newtonian equation md 2 r/dt 2 = F. In that 
time quantum mechanics has not been found. Whether the Newtonian equation can keep unchanged under 
time reversal depends on the form of force F. Only when F is conservative, it dose. In general situations 
when F is relative to time t and momentum p, the Newtonian equation can't keep unchanged under time 
reversal in general. In the common statistical systems composed of charged micro-particles, the interaction 
forces between micro-particles are electromagnetic forces. In the systems composed of neutral atoms and 
molecules, atoms and molecules can be regarded as electromagnetic dipole moments and quadrupole moments 
for their deformations caused by interactions and interaction between them can considered as one between 
electromagnetic dipole moments and quadrupole moments. In the paper, we only discuss interaction between 
charged particles, but the principle is the same for neutral atoms and molecules. 

If the retarded interaction is not considered, the Lorentz force between two particles with charges q and 
q' as well as velocities v and v' is 

c? iq'r , qq'v' x (v x r) 

— + ^2 W 

In the current statistical physics, we use the formula to describe interactions between micro-particles. Because 
the force is unchanged when v — > — v and v' — > —v 1 under time reversal, it is impossible for us to introduce the 
irreversibility of time reversal into the theory to solve the problems in non-equivalent statistical mechanics. 
However, it should be noted that according to special relativity, instantaneous interaction does not exist. 
Electromagnetic interaction propagates in the speed of light. So retarded interaction should be considered 
between micro-particles in statistical systems. It is proved below that after retarded interaction is considered, 
the Lorentz forces can not keep unchanged under time reversal again and a rational dynamic foundation can 
be established for classical statistical mechanics and the reversibility paradox can resolved well. 

Let t',r',v' and a' represent retarded time, coordinate, velocity and acceleration, t, r, v and represent non- 
retarded time, coordinate, velocity and acceleration. A particle with charge qj, velocity v'j and acceleration 
a'j at space point rj(t') and time t' would cause retarded potentials as follows at space point fi(t) and time t 

v = — $_ — 4- - — q -&, — (2) 
(i - c(i - -v 1 )^. 

In the formula f^(M') = n(t) - rj-(t'), r'^ =| ^ |, n'^ = ^/r'^. Let v' jn = v] ■ n' jn , a' jn = a'j ■ n'^ the 
intensities of electromagnetic fields caused by j particle at space point u and time t are 

e ad-^-fi) ^x [( ^ fi)x gj | 

(1 - -^fr 12 c 2 (l - -^fr 1 ^ 

We call the forces as the retarded Lorentz forces in which retarded interaction has been considered. Suppose 
there are N particles in the system. Suppose the i particle with charge qi at space point r,(t) at time t 
moves in velocity i>i, the retarded Lorentz force acted on the i particle caused by the j particle is 

fi _ yy2 & c c 

K c > ij 

c 2 (l-^)%- 

-^(l-^)-a'(l-^)(l-^)} (4) 
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In the formula, vi n — - Vi, v'j n — ftij - Vj, a'j n = n'^ -a!y In order to do rational approximate calculation, let's 
accumulate the magnitude order of acceleration. According to formula a = v 2 /r', a/c 2 has the magnitude 
order of \jr' . Suppose hydrogen atom can be regarded as a harmonic oscillator with amplitude b and 
angle frequency oj, so the acceleration of oscillator is bu; 2 . The magnitude order of energy of hydrogen 



atom is E — hc/X — mv 2 — mb 2 uj 2 . Suppose the wavelength of photons emitted by hydrogen atom is 
A = 4 x 1CT 7 M, it can be calculated that a/c 2 = 1/r' ~ 60. But the magnitude order of distance between 
atoms is r ~ 10 _10 M, 1/r ~ 10 10 . According to Eq.(4), the first item of the retarded Lorentz force directs 
ratio to 1/r 2 . So for the neighbor interaction, we have a/c 2 r = 1/rr' = 6 x 10~ 9 /r 2 << 1/r 2 , the item 
containing acceleration can be omitted. For the distance interaction, a/c 2 > 1/r 2 , the items containing 
acceleration can not be omitted On the other hand, according to the Maxwell distribution of velocity, the 
average speed of hydrogen atoms is v — yjBkT firm. Under common temperature T = 300K, we have 
v = 6.3 x 10 6 M, v/c = 2.1 x 1CF 2 . So the item containing v/c is much bigger than the item containing 
acceleration, even the items containing v 2 /c 2 are bigger than that containing acceleration in near neighbor 
interaction. So we should retain the items containing v/c, v 2 /c 2 and va/c 3 , but omitted high order items, 
then write Eq.(4) as 

Fru = "^M 1 + — + — -^t- — ) - -(! - — + — )1 
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For the convenience of discussion later, we write F' Ri - = F^- +F-p F^ - = qiqjn'^/r'- 2 represent conservative 

part and F[- represents non-conservative part. The total force acted on the i-particle caused by the other 

particles is F' Ri — X^-^oij + S-^/j = ^oi + ^ or ^ e macro-systems composed of large numbers of 
neutral atoms and molecular, we can consider atoms and molecular as electromagnetic dipole moments and 
quadrupoles and obtain the retarded Lorents forces in the same way. 

On the other hand, according to classical electromagnetic theory, there exists the radiation damping 
forces G acted on the accelerated particles. The radiation damping forces are relative to the acceleration of 
acceleration. Suppose I is the dimension of the particles, particle's speed v << c, acceleration a << c 2 /I, the 
acceleration of acceleration d << c 3 // 2 , when the distribution of particle's charge is with spherical symmetry, 
we have G = net, k — 2q 2 /3c 3 . 

In classical electromagnetic theory, the Hamiltonian equations of motions can be described by using the 
space coordinates and common momentums Pi, or by using the spaces coordinates and canonical momentums 
p z i with relation p Z i = pi + qiAi/c. Both are equivalent. In statistical mechanics, it is more convenient to 
use common momentums. The accelerations and the accelerations of accelerations of particles should be 
regarded as the function of coordinate, momentum and time. On the other hand, the acceleration of the 
-particle is caused by other particle's interactions, and it is relative to the retarded speeds and distances, i,e., 
al i — a' i (fi,Pi,r'j,p'j), a i = a i {r i ,pi,r 1 ^'^^). So the total Hamiltonian of the non-conservative system composed 
of N charged particles can be written as H = H + H' , in which H is the conservative Hamiltonian 

i 1 i<j 

Here C/ij(r^) is the conservative interaction energy, H' is the non-conservative Hamiltonian 

H' = ^ Ui<?i,t) + £ UiMvPiJi) (7) 

i i<j 

Ui(fi, t) is the interaction caused by external force, Uij (r ^ , Pi , p*j) is the interaction caused by non-conservative 
force. The motion equation of the i-particle is 



pi TT pi TT 



Q' ia = F eirT (fi , t) + F[ a (n ,fj,Pi,ifj) + G' ia (n ,fj,Pi, p 3 ) (9) 

In the formula, is conservative force Q' ia is total of non-conservation force, F ei(7 is external force, F( a is 
total retarded non-conservation force, G' ia is radiation damping forces. Because acceleration of acceleration 
of a particle is caused by other particles, we can write G' ia = J2j^i ^'ija m general, G' i j rT contains the j 
-particle's inference on the acceleration of i -particle's particle. We can express f, r', v' and a' by r, r, v and 
a as shown in Eq.(35)-(38). Because when t — > —t, we have v — > —v, a — * a, so that the retarded Lorentz 



forces and the radiation damping forces can not keep unchanged under time reversal, and irreversibility 
would appear in theory. These will further be shown in detail later. 

Two basic hypotheses can be regarded as the foundation of classical statistical physics as shown below. 

1 . The interaction forces between charged micro-particles in the macro-systems are the retarded Lorentz 
forces and the radiation damping forces showing in Eq.(9). For neutral atoms and molecular, the interac- 
tions forces can be considered as the retarded Lorentz forces between electromagnetic dipole moments and 
quadrupoles. 

2. A classical statistical system can be described by the normalized distribution function of ensemble 
probability density p = p{x ia ,p ia , t). The average value of the physical quantity u in the ensemble is 

N 3 

(10) 



p la , t)p(x ia , p ia , t) Y[ Q dx la dp 

-> i=l er=l 



We will establish the dynamic equations of classical statistical physics based on these two hypotheses, 
then discuss the problems of equivalent and non-equivalent states below. 

3. The basic dynamic equation of classical statistical physics 

After the retarded Lorentz force and radiation damping force are introduced, the time rate of change of 
the distribution function of ensemble probability density p is 

1(7 1(7 

By using Eq.(6) and the continuity equation 

1(7 

Eq.(ll) can be written as 

1(7 1(7 

From Eq.(4) we can obtain 



OF' 



dpt, 



= (14) 



Put Eq.(14) into Eq.(12), we get 



dp_dp s^ p^dp , , dp i _ n \- dG l 

Tt at + ^ t dx~ + {Fei ° + Foi ° + Fi ° + G ™ ] dp- ] -~ p ^^p~ (15) 

1(7 1(7 

But dG' ia /dri a ^ in general, for example, for a charged oscillating dipole, we have x = Asin(ujt + 5), 
p = rax = Amuo cos(uot + 5), a = —Aiv 3 cos(u;t + 5) = —uj 2 p/m, so we have dG/dp = —nuj 2 /m ^ when the 
retarded effect is not considered. Eq.(15) can also be written as 

at + ^ [ mdx- + {F - + + F - ] dp- ] - -^p— (16) 

1(7 

The Eq.(15) or (16) are just the motion equations the distribution function of ensemble probability 
density p satisfy after the retarded Lorentz forces and radiation damping forces are introduced, it can be 
regarded as the basic equation of the classical statistical mechanics. When F[ a = G' ia = 0, the equation 
becomes the current Liuve equation. Because F' ia and G' ia can not keep unchanged under time reversal, 
Eq.(15) and (16) also can't be unchanged under time reversal. 



4. The statistical distribution of equivalent states 



Let's first discuss the definition of equilibrium states and the hypothesis of equal probability in statistical 
physics. The equilibrium states in thermodynamics have strict definition. The equilibrium states are the 
states that the system's macro-natures do not change with time without external influence. But in classical 
statistical theory the equilibrium states have no independent and strict definition. For the equilibrium states 
of isolated systems, the definition depends on the equal probability hypothesis. The hypothesis is that for 
the equilibrium states of isolated systems, the ensemble probability density is a constant between the curved 
surfaces of energies E and E + AE in phase space. So in current theory, the equilibrium states of isolated 
systems are equivalent to the states of equal probability. For the equilibrium states of non- isolated systems, 
there is no strict definition now. Some documents consider them as the states that the average values of 
physical quantities do not change with time. But this kind of definition depends on what kinds of physical 
quantities are used. So it is improper, for the equilibrium is only the nature of system itself, it has nothing 
to do with other physical quantities. As for the hypothesis of equal probability, it is a priori hypothesis 
and can't be proved in theory. The reason for the hypothesis to exist may be based on so-called " the law 
of non-difference". The hypothesis had got much criticism since it was put forward. In order to provide 
the hypothesis a rational base, Boltzmann raised the ergodic theory, proving that as long as a system was 
ergodic, the hypothesis of equal probability would be tenable. However, the study shows that the evolutions 
of systems can't be ergodic generally ly2 \ So the ergodic hypothesis can't be used as the foundation of 
classical statistical physics and we have to look for other way out. 

For this purpose, let's discuss the definition of equilibrium states in statistical physics at first. We 
define the equilibrium states of isolated systems as the states that there arc no external forces acting on 
the systems and the ensemble probabilities do not change with time. This definition coincides with the 
definition of thermodynamics of equilibrium states. According to the definition, for isolated systems we have 
F e j(r, t) = 0. When the equilibrium states are reached, we have dp/dt = and dp/dt — 0. So in the light of 
Eq.(15) in equilibrium states, we have 



They are just the equations that the ensemble probability density functions satisfy in the equilibrium states. 
There are many kinds of equilibrium states of isolated systems as shown below. 

1. The particles are free in the systems with = F- = G\ = and the ensemble probability density 
satisfies dp/dx^ = and dp/dpi a = 0. In this case, the ensemble probability density function p = constant, 
i.e., p has nothing to do with x irT , p ia and t. This is just the state of equal probability or the so-called 
micro-canonical ensemble in the current theory. 

2. Particles in systems are only acted by conservative forces, F^ ^ 0, F[ = G\ = 0, and dp/dxi a ^ 
dp/dpia 7^ 0. The retarded effect is not considered. In this case, we have 

dp p ia , dp_ p = _dp_dHo _ _dp_dHa_ = Q 



dx ia rrii dpia dx irT dp ia dp ia dxi 

In this kinds of equilibrium states, the ensemble probability density functions do not change with time, but 
the probabilities are different in the different points of phase space, or the equivalent states are with unequal 
probabilities. It is easy to be verified that the solution of Eq.(18) is 

p = Ae-< m = e~^ E (19) 

Here A, (3 and ip are constants. E = Hq is the total energy of the system. In fact, according to Eq.(19), we 
have 

dp = _p dHo_ dp_ = _p dHo_ 

dx ia dx ia dpi a dpi a 

Put it into Eq.(18), we can prove that it is the solution of Eq.(18). Eq.(19) is just the canonical ensemble in 
the current statistical theory. It is known that this result has nothing to do with whether or not the retarded 
Lorentz forces and the damping radiation forces are introduced. Based on the current Liouville's equation 
and the definition of equilibrium states of isolated systems, we can get it directly. But unfortunately, this 
simple method is neglected by the current theory. Perhaps it is just this neglect, the hypothesis of equal 
probability became necessary and possible to exist. 

If the particles in the system are identical, and the force acted on each particle only relative to its 
coordinate, or the particles are independent each other without interactions between them, we have E = 



Ei, Ei = Ti + Ui. Here Tj is particle's kinetic energy, Ui is potential energy, Eq.(19) becomes 



p = Ae-^ El+ - Ei+ - +E ^ = Ae~ 0El ■ ■ ■ e~ 0Ei ■ ■ • e^ E " (21) 

The formula is just the distribution of near-independent subsystem. By the integral over while space and 
the normalization of the p function, we get 

/+oo r+oo r+oo r+oo 

pdil = / A^e-^dn, ■ ■■ A^e-^dSU ■ ■■ A l ' N e"^ dCl N = 1 (22) 
-oo J — oo J — oo J — oo 

Because the particles are identical, the forms of their energies E{ are the same, so Eq.(22) can be written as 

/+oc z' + oo p-\-oo 

pdn = ( A 1 ' N e-^ Ei dVL i ) N = 1 or / A^e^dQi = 1 (23) 

-oo J — oo J —oo 

Let A x / N = e~ a , the formula can be written as the form of sum 

fi^i = E e" a - /3£i AO i - 1 fiASli = e- a -^A^ (24) 

i i 

This is just the distribution law of the Maxwcll-Boltzmann. or the distribution law of most probable value. 

If the particles are free without potential forces, we have Ui = and Ei — > mv 2 /2. So the formula above 
becomes the Maxwell distribution law of velocities. 

mv 2 

fdQ. = Be~ii<T dxdydzdv x dv y dv z (25) 

Here KT = 1/(3, B = e~ a . In this case, the distribution function has nothing to do with space coordinates, 
showing that the distribution of particles in space is uniform. 

When all particles have the same speed in system, according to Eq.(25), we have / = constant and the 
system becomes micro-canonical ensemble. So micro-canonical ensemble represents the systems composed 
of free particles, in which particle's speeds are the same and the space distribution is even. This kind of 
systems, of course, has no significance in practice. 

On other hand, if there are two systems with the numbers of particles N\ and N2, energies E\ and E2 
individually, we can also get the grand canonical ensemble as done in the current theory. 

So it is obvious that we can get all kinds of equivalent distributions from canonical ensemble, including 
micro-canonical ensemble. But we can't obtain canonical ensemble from micro-canonical ensemble. The 
situation is just opposite comparing with traditional theory According to the current theory, canonical 
ensemble is deduced from micro-canonical ensemble, representing the equilibrium states of small systems 
to contact with big heat sources. So canonical ensemble describes the equilibrium states of non-isolated 
systems according to traditional statistical theory. But according to the Liouville's equation, canonical 
ensemble describes the equilibrium states of isolated systems. So micro-canonical ensemble contradicts the 
result of the Liouville's equation. On the other hand, micro-canonical ensemble is only one of equilibrium 
states. The other equilibrium states are not with equal probabilities. So it is improper to consider micro- 
canonical ensemble as the foundation of equilibrium state theory. For the systems only acted by conservative 
forces without considering retarded effects, the canonical ensemble is more foundational. For the general 
statistical systems, the revised Liouville's equation is fundamental. It can be seen from discussions above 
that it is more simple and rational to discuss the problems of equilibrium states from the Liouville's equation 
directly, no matter on physical concepts or mathematical calculation. 

3. When the non-conservative forces are considered and the system satisfies 



dpP -+^{F^ + F la )] = , (26) 



dx ia rrii dpi 



This is other equilibrium state. There are other equilibrium states. 

For non-isolated systems, we can definite the equilibrium states as the states that the forces acted on 
the systems and the ensemble probabilities functions do not change with time. For example, the equilibrium 
state of a system located in the gravitational field of earth. The distribution function can still be described by 
Eq.(19), but in this case total energy E contains potential energy. In all theses equilibrium states, ensemble 
probability densities change with the coordinates of phase space, but do not change with time. 



Therefore, it is clear that we can deal with all the problems of statistical physics well from the revised 
Liouville's equation directly. It is unnecessary for us to introduce some extra hypotheses such as the equal 
probability hypothesis. At present, chaos theory is also involved when we discuss the problems of foundation 
of statistical mechanics. Based on conservative interaction, thought chaos theory may be relative to the ori- 
gin of random nature of statistical system or the hypothesis of ensemble, it has no effect on the foundational 
form of motion equation of statistical mechanics and has nothing to do with the origin of irreversibility in 
macro-systems. So it is impossible to us to solve the fundamental problems of statistical physics depending 
on chaos theory. In fact, it is enough for us in principle to solve all problems in statistical physics based on 
the revised Liouville's equation. By the method provided in the paper, the united description of equilibrium 
states and non-equilibrium states without using the concepts of ergoden, coarse graining and mixing currents 
and so on, which are easy to cause dispute. 

5. The statistical distributions of non-equivalent states and BBGKY series equations 
For non-equivalent states, by the direct integral of Eq.(15), we have 

,(*) = ex P [- £ / d -9k dt] = ex P [- E E E / ^M^^iM^M^] (27) 

i(T a—1 i—1 j^l 

In order to let the integral is possible, the retarded quantities should be expressed by the non-retarded 
quantities or opposite. It is known that the relation between retarted time and particle's distance is 



ra = c(t t 1 ) = yJ[ Xi (t) x' 3 {t')f + [ Vl (t) y'(t>W + Wt) z' 3 {t'W (28) 

Suppose the function relations Xi = x,(t) and x^ — x'j(t') are known, it can be obtained from the formula 
above 

t = f[x' j (t'),y' j (t'),z , j (t'),t'} or t 1 = f[x i {t),y i {t),z i {t),t] (29) 
From the relations, we have 

dt' dx' 3 dt' dy'j dt' dz^ dt' y dt 1 3 ' v ' 

Put it into Eq.(27), we get the probability function expressing by retarded time 

cr=l i=l j^i la 

On the other hand, by using relations x, = Xi(t), yi = yi(t), Zi = Zi(t) and Eq.(29), we have 

r'^ (t, t 1 ) = c[t - f'{xi, y u z u t)} = r'^ (t) 

& j =if j [f'{xi,y i ,z i ,t)] = if j (t) a' j =a' j [f'(x i ,y i ,z i ,t)]=a' j {t) (32) 

Put them into Eq.(27), we get the probability function expressing by non-retarded time 



(33) 



When the particle's speed v « c, the retarded distance r^(i') in the retarded time t' can be replaced 
approximately by non-retarded distance rij(t), i.e., we can let 

t' =t-r' j (t')/c^t-r ij (t)/c (34) 

r'^t, t>) = | m ?. (f ) | = | m -r^[t- r[ 3 {t')/c] |-| n(t) - - r«(t)/c] |= r'^t) (35) 

It is noted that r[j(i) ^ fij(t), for r^-(t) is the approximate retarded distance, but Tij(t) is not the retarded 
distance. In this case, we can develop retarded quantities into series in light of small quantity nj/c. By 
relation vj = dfj/dt, we get 

#^(t,f) ~ n{t) ?.{t r ijlc ) = nit) f 3 (t) + r -^-vj(t) r ^§-d, + r ^k 3 + ■ ■■ 



- n 3 (t) + «j (t) - -^-oj + -^3- a i + • • • ( 36 ) 
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dj(f ) ~ ^(t) - -^-(t) + + • • • S$(f ) ~ a,(i) - -^-(t) + • • • (38) 



From Eq.(36) and (37), it is obvious that rjy(t') ^ fji(*')> r 'i 3 (t') r ji(0- The relations below are useful in 
later calculations. 

J_ ^ ±{i - v JlL + (i - ^ )ri . _ (i _ ^H) r 2 + . . .} (39 ) 
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The BBGKY series equations are discussed now. The volume element of phase space of the i-particle is 
written as dCli = U<j=idxi a dpi a and the reduced distribution functions of ensemble probability densities are 

fs(xia,Pla ■ ■ ■ x sa ,p s<T ,t) = J ' p(xi a ,pi a ■ ■ ■ x nrT ,p na ,t)dQ. s +i ■ ■ ■ dfl N (42) 
In the current theory, by considering the nature of identical particles, the BBGKY series equations are 

r\ p 3 5 <-\ p i~\ p 3 5 5 

ot i r-f m dx ia dp l<y ^ ^ dp lcr 

a—Y i—i a—I i=l j^ti 



(- V - S ' : EE f Fis+ia^dSl^ (48) 

<T=1 i=l J 



The equations are equivalent to the Liouvillc's equation. Considering the fact that F eia only relative to the 
i -particle's coordinates but Fg i(T , F' ia and are relative to both the i and j particle's coordinates after 
the retarded Lorentz forces and radiation damping forces, we can get new BBGKY series equations in the 
same way 

r\ /> 35 C\ £ €S £ 355 O-f 

1 + E Eelfc + ^) + E E E^ + ^ ' 



+(^-5)11/(^+^)^+1 

<r=l i=l J ° Pla 



(44) 



5=1 !=1 j^i <T=1 i = l • 



In the formula, FQ iS+lt7 , F! s+1 a and G' iS+1 a represents the j = S + 1 items. Because iV — 1 ~ N, the first 
equation with S — 1 is 

= -N J F[ 2 ■ V p ~J 2 d 3 r 2 d 3 p 2 -N J V fl ■ (G' ij f 2 )d 3 f 2 d :i p 2 (45) 



In the formula, Fg 12 is the conservative retarded force and F[ 2 is the non-conservative retarded force. By 
using Eq.(36) — -(41), F\ 2 tcan be expressed by the non-retarded quantities as 

Foa = ^ = ^ + K l3 = -V* • U( nj ) + K l3 (46) 

^ - —3—^1 + "Til- + -Z=— K ij3 + — — K ij* ( 47 ) 

Ktfi = f- + + "^-C 1 - -7-) r « - ^ - -7-)^ + • ' ' ( 48 ) 
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(49) 



It is obvious that after the retarded quantities are expressed by the non-retarded quantities, the conservative 
forces would become the non-conservative forces. In the same way, G' io can also be expressed by the non- 
retarded quantities. By using Eq.(46), Eq.(45) can be written as 



it + m ' V?1 h + / ' V?1 h ~ N I V?1 U{ri2) ■ Vpl f "' :i r 



J VfiI7(ri2) 



= -N J (K 12 + F{ 2 ) ■ WpJ 2 d 3 f 2 d 3 p 2 -N J V ft • (G' 12 f 2 )d 3 r 2 d 3 p 2 (51) 

The left side of the equation is the result by the Liouville's equation. The right side is new by the paper's 
revision. Similarly, the second equation of the BBGKY series with S = 2 can be written as 

+ - ■ V*/2 + - • V r - 2 / 2 + F el • V ft / 2 + F e2 • V p - 2 / 2 - VfiE^na) ■ V^/ 2 
at m m ri 

-V r - 2 (7(r 12 ) • V p - 2 / 2 - N j [V Pl U(r 13 ) ■ V ft / 3 + V ?2 (7(r 23 ) • V p - 2 / 3 ]d 3 r 3 d 3 i? 3 

= -V fi • (G' 12 / 2 ) - V p - 2 • (G' 2 J 2 ) + K 12 • + £ 21 • V p - 2 / 2 

-TV J [{K 13 + F{ 3 ) ■ V Pl f 3 + (K 23 + F 23 ) ■ VpJ 3 ]d 3 f 3 d 3 p 3 

-N j [V ft (G' 13 / 3 ) + V p - 2 (G 23 .f 3 )]d 3 f 3 dV 3 (52) 

The definition of Kij can be seen in Eq.(47). We will discuss the motion equations of hydromechanics me- 
chanics by using them below. 

6. The motion equations of hydromechanics mechanics 

How to reduce the motion equations of hydromechanics from statistical mechanics is still a unsolved 
problem now. Because the current statistical mechanics is based on the Liouville's equation, it is only 
suitable for the equivalent states of conservative systems or the ideal fluids. The dissipative phenomena of 
real fluids just as heat conduction and viscosity and so on can only be explained rationally by introducing 
non-conservative forces. Let's discuss this problem now. By writing the right items of Eq.(50) into the 
form of hydromechanics and adding them into the current results, we can get the results after the retarded 
Lorentz forces and radiation damping forces are introduced. According to the current theory, we definite 
the normalization of functions f\(f\.pi,t) and f 2 {f\.p\,r 2 ,p 2 ,t) as 



V 



= f MruPi^nd 3 ^ (53) 



Vo 2 = I h (n ,Pi,r 2 ,p 2 , i)d 3 fid 3 pid 3 f 2 c?V2 



(54) 



In the formula po is macro-mass density, V is velocity of fluid, Uk are dynamic energy density, u v is potential 
energy density 

Po{ri,t) = —J /i(fi,pi,i)d 3 pi V(n,t) = — —J pi/i(ri,j?i,i)d 3 pi (55) 

= / (P- my) 2 /i(ri,pi,i)d 3 pi (56) 

- -^(^) 2 / U(f 1 ,r^)f 2 (n,P2,f' 2 ^ 2 ,t)d 3 p 1 d 3 r 2 d 3 p 2 (57) 



2po Vo' 

The equilibrium equation of mass density or the continuity equation is discussed firstly. Producing 
Eq.(50) is produced by mN/Vo, then the integral about d 3 p\ is carried out The right side of the equation 
can be written as 

N 2 /\.t - _ . 



PoV 



r _ N 2 f ^ 

J (K 12 + F[ 2 ) ■ VpJ 2 d 3 p 1 d 3 f 2 d 3 p 2 --^y o J V fi • (G' 12 / 2 )dVid 3 r 2 dV 2 (58) 



According to Eq.(14), we have • F{ 2 = 0. Because i4Ti 2 has nothing to do with pf, we also have 
Vpj • i?i2 = 0. By using the boundary conditions of probability density functions, Eq.(50) is zero. So the 
continuity equation of hydromechanics is the same as the current form with 

^+V-p Q V = (59) 

In order to get the motion equation of hydromechanics, we product Eq.(50) by Npx/Vo and take the 
integral about d 3 pi, then get from left side of Eq.(50) 

Po(^ + V-VV)+V-(Pk + K)=PoTi (60) 

In the formula, Ti , Pk and P v are the average external force, kinetic energy and potential energy tensors of 
unit masses individually ( 5 ) 

N t s - N r * - m ^)(pi - ^ 



PoV J V J m 



ft = - 2 ( V YJ o dXj ^^^i 2 d^d 3 p 2 (62) 



For the right side of Eq.(50), because • [K\ 2 + F[ 2 ) = 0, we get 

N 2 f N f 

-— j Pi(K 12 + F[ 2 ) ■ Vpj2 d 3 Pid 3 r2d 3 P2 -yj ftV ft • (G' 12 / 2 )dVid 3 r 2 dV 2 

N 2 f N 2 f 

= y£j M K ™ + ' VpJi^Pid^fc +yrl hG\ 2 • Vp l pid 3 p 1 d 3 r 2 d 3 p2 

N 2 f -> N 2 I' 

= ~rj~ / /2(^i2 + F\2 ) d 3 Pi d 3 f2 d 3 p2 + — / f2G' 12 d 3 p 1 d 3 r 2 d 3 p 2 (63) 
"o J Vo J 

On the other hand, by using Eq.(36) — (41), Eq.(5) can be expressed by the non-retarded quantities 

ft, <jj(ljTjj ( Mf u ij n . QiQj^ij n ( Mj CL ij ri in A \ 

*H ~ —3—^1 + rr 2 + o, • Qij3 + o Qijl ( b4 ) 

* ij Cl ij L '%J c 

Let Rijk = Kijk + Qijk, we can write K\ 2 + F[ 2 as 

# p, _ gi 92^12 D gig 2 ^2 D gig 2 a 2 p gig2«2 D , . 

^12 + ^12 — 5 -n-121 H o ^122 H 5 ^123 H 3 J1124 (,OOj 

rj 2 crj 2 c^r 12 e 3 



In the formula, Rij k does not contain retarded quantities again. By the relation qiq 2 ff 2 = — Vf 1 J7(ri2), we 
have 

/V 2 (■ _ 

-yr J (K 12 + Fl 2 )f 2 cPp l cP?2<f i & 

= TT [ ^P 1 U{r l2 )R l21 f 2 d 3 Pid?r 2 d 3 p 2 -^- / ^^fli22/2dVid 3 r 2 d 3 p 2 
Vb 7 Vb y cr 2 2 

N 2 f qiq 2 a 2 , ^3 - ,,3 - j3 - ^ 2 f <Hl2a 2 3 ^ 3 ^ 3 ^ . . 

-77- / — o— #123/2^ Pid r 2 <rp 2 - 77- / — 2—Ri24f2<rpi<rr2<rp 2 (66) 
Vb 7 crf 2 Vb y cr 2 2 

By using formula ( 5 ) 

Jf(n,f2)V ri U(ri2)d 3 f2 = y ri -{-lj\\ J dM' 2 ^ ^g 2) /(rl + (1 - AK 2 , n - Af/ 2 )} (67) 

/ V fl U{r 12 )Ri2if2d s pid s r2d 3 p2 = -V* • P s (68) 

^ = ^ / lrfA / ^5 a ^ ( ? a) ^i[^ + (i-A)^,ft,n-Ai? a ,g,,t] 

^OJo 7 r 12 « r 12 

x/ 2 [ri + (1 - \y{ 2 ,Pi,n - K 2 ,P2, i]rfV/ 2 rfVirfV2 (69) 
We call as dissipative energy tensor. Let 

f 2 = f ( ^ij 122 + l^ Rl23 + %R 124 ) h d 3 p ld 3 r 2 d 3 p 2 (70) 

cpoVb y rf 2 cri2 c 2 

represent the average dissipative force acted on unite mass and 

n 2 r 

r 3 = —tt / f 2 G\2d Z Pid z r 2 d z p 2 (71) 



we have 

TV 2 

r2 /■! 



PoVb 

represent the average radiation damping force of unite mass, the motion equation of hydromechanics is 

Po{ % + V-VV) + V-(P k + P v + P s ) = p (f i + f 2 + f 3 ) (72) 

Producing Eq.(50) by N(p\ — mV) 2 /2mVo and taking the integral about d 3 pi, we can get the equivalent 
equation of kinetic energy. Form the left side of the equation, we have 

+ V • ( Po u k V + J k ) = -P k : VV (73) 

Here J k = J kx + j k2 with 

N [p-mVfa -mV) 2 . j3 ^ 
Vb y m 2m 

J k2 = f dX ^ 2 ' (P V dU } r 'I 2) h{ri + (1 - A)*? 9 ,ft,^ - Xr^Md^Pid^ (75) 

JO ™12 « r 12 

Form the left side of the equation, we get 

n 2 r - - 

~2mV J ifl ~ mt/)2V w ' [{K + F l2)h]d 3 Pid 3 f 2 d 3 p 2 
N 



J (p! - mV ) 2 V P - • (G' 12 / 2 )dVirf 3 ^2rfV2 

n 2 r 



2mVb 

r2 



+ I h&12 ' {Pl ~ m V) d3 Pi d3 r2d 3 P2 (76) 
By using Eq.(65), the first item of the formula above can be written as 

/ / 2 (¥^m + 4-^22 + -#-^123 + ^i?i 24 ) • (pi mV)d?nd?r 2 d*p 2 (77) 

J r \2 cr 12 C r 12 c 

By using Eq.(67) again, the first item of the formula above can be written as 

v • J k3 = 4r / hRui^^P 1 ■ (pi - mv^p^d 3 ^ 

mVo J r? 2 

? iV 2 f 1 „ f ■ (p x - mV) dU(r'{ 2 ) „ ^ , xW/ ^ ^ ^ , 



12 12 



xf 2 [f 1 + (1 - AK 2 ,Pi,n - \r{ 2 ,p 2 , t^l^p^pz (78) 

Jfc3 can be called as the dissipative fluid of kinetic energy of unite mass. The second, third and forth items 
in Eq.(77) and the last item in Eq.(76) can be put together and written as 

Po<?k = / / 2 (^Qi22 + — 0i23 + ^Qi24 + — G' 12 ) ■ (n - mV^p^fzd 3 ^ (79) 

cmVo J r{ 2 cr 12 c 2 q x q 2 

We can call <7fc as the dissipative kinetic energy generation of unite mass. In the current theory, there is no 
this item. It appears only when the retarded Lorentz forces and damping radiation forces are introduced. 
Let Jfe = Jki + Jk2 + Jk3, thus Eq (73) can be written as 

^^- + V-(p u k V + Jk) = -3k ■ VV + p a k (80) 

Producing Eq.(51) by N 2 U(ri2)/V 2 and taking the integral about d 3 pid 3 r 2 d 3 p2, the equivalent equation 
of potential energy can be obtained. Form the left side of the equation, we can get 

+ V • (p u v V + J v ) = -I : W (81) 

Here J v — J v \ + J V 2 with 



Jul = \{^f I ^-^U(r ij )f 2 d 3 p 1 d 3 f 2 d 3 p2 (82) 

r, 2 = - i rfA ^% {p l -p 2) du } r J 2) f2(n + (i - xy; 2 ,p 2 ,n - a^^MWi^ m 



Jo 



Amr'l 2 dr 12 
The right side of the equation is 

T2 



n 2 r n 2 f 

y5 J U(r 12 )K 12 ■ Vp 1 f2d 3 pid 3 r 2 d 3 p 2 + j U(r 12 )K 21 • V p - 2 / 2 d 3 pirfV 2 rfV 2 

n 3 r - - 

-^5- J U(r 12 )[(K 13 + F{ 3 ) ■ V P J 3 + (K 23 + F^) ■ VfJ^prd 3 ^ 3 ^ 3 ^ 3 ^ 

N 2 f - N 2 f 

~y2 J f / (^i2)V P - 1 • (G' 12 / 2 )dVid 3 ^rf'V 2 -ygj U(r 12 )V A ■ (G 21 / 2 )rfVid 3 r 2 d 3 p2 

-JV J (7(r 12 )[V p - • (G' 13 h) + V p - 2 • (GQf 3 )]<Pp\ < P?2<Pp2<P?3<P& i (84) 

As mentioned before, we have Vp t ■ = 0, • K\- = 0. By considering the boundary condition, all 
integrals are zero. So the potential energy items have no reversion after retarded Lorentz forces and the 
radiation damping forces are introduced. 

We have to consider the contributions of non-conservative dissipative energies. It is known in electro- 
magnetic theory that when the common coordinates and momentums are used to describe particle's motion, 



interaction energy does not contain magnetic field's action. The retarded interaction between the i -particle 
and the j -particle is 

ulMirt) = , n qiq \ n = ^(i + — + 4) (85) 



U[j is asymmetry about the indexes i and j after the retarded effect is considered, so the non-conservative 
and retarded total interaction can be written as 



JV JV / , 2 N N 

w^) = \ E E f-tf + ^) = ^ E E 7^ 

i=i jyi tj i=1 j-^! y 



x{ ^_^ (1 _^ ) + ^ (1+ 4^ ) _3^-y tJ+ 2u,.^ } 
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2c 3 v 3c ' 2c 3 3c 4 



Producing the second equation of the BBGKY series equation by N 2 U'(r' 12 ,pi,p' 2 )/2V Q 2 an d taking the 
integral about d 3 p\d 3 r 2 d 3 p 2) we get 

8 N 2 f N 2 f 

dtZV 2 j V '^^ d ^ 2d ^ 2 + J U '& ■ V ^ +P2 • ^r 2 ).f2d?p l d 3 r 2 d 3 p 2 

N 2 f 

2 /■ 

^ y U'[V Fl U{r 12 ) ■ V ft + V ?2 [/(r 12 ) • V p - 2 ]/ 2 d 3 pid 3 r 2 d 3 i? 2 
^ y f/'[V ?1 [/(r 13 ) • V fi + V r - 2 C/(r 23 ) • V^]hd 3 ^ 3 r 2 d 3 p 2 d S r 5 d 3 p 5 

n 2 r 

_/ ^'[Vp. • (G' 12 / 2 ) + V p - 2 • (G 21 / 2 )]rfVirf 3 r 2 rf 3 p 2 
N 2 f 

+ 2V 2 J U ' {Kl2 ' V?1 + K21 ' V P2)f2d 3 pid 3 f 2 d 3 p 2 
y C/'[(^i3 + ^'3) • Vp-Ja + (#23 + F^) ■ V p - 2 / 3 ]dVirf 3 r 2 dV 2 d 3 r 3 dV3 

TV 3 r 

'W 2 j U ' [Vpl ' (G ' 13/3) + Vft ' ( G 23/3)]rfVirf 3 r2dV2d 3 r 3 rf 3 P3 (87) 
The calculation is verbose and we only give out the result of the first item in the right side of Eq.(86) 
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~2V r 
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~2Vr 



TJ q\q2 ( v 2n . vi n q\q2,r\2-v Y f 21 ■ Vi, qiq 2 r 12 ■ {p 2 - Pi) 

Ui ^ { ^- ~^r ] = 2^ — 15 — (88) 



Similar to Eq.(56), we definite 

IN 2 , , 

u s = -z — ttj / Uihd^d^d 3 ^ (89) 
zp v 



Jsl = ^ l^^Vlhd 3 pid Z ?2d 3 p 2 



as the energy density and energy fluid density of non-conservative interactions of unite mass and have 

Pou s V = i(^) 2 J VUxhd Z Pid 3 ? 2 d 3 p 2 (91) 

V r - • (p u s V + Xsi) = V* J ^U!f 2 d 3 p!d 3 r 2 d 3 p 2 } 

= It^ /Vi- ' V n/2 + h~ ■ V ?1 C/i)dVirf 3 r 2 d 3 p 2 (92) 
2 Vq J m m 



By considering the boundary condition, we can get 

N 2 f N 2 { 

j U\Vi ■ Vr 2 .f2d 3 pid 3 r 2 d 3 p 2 = -^72 J hP2 ■ Vf 2 Uid 3 pid 3 f 2 d 3 p 2 (93) 



2mV r 
Therefore, we have 



N 2 f 

^—^2 J Uitfi ■ V ?1 +P2 ■ Vf 2 )f 2 d 3 p 1 d 3 r 2 d 3 p 2 

V* • (poUsV + J sl ) - I f h ■ V ft Uid 3 pid 3 f 2 d 3 p 2 (94) 
^ > / o J m 



For Eq.(88), because 

d {xi - x 2 ){p 2x - pi x ) = p 2 x - Pis _ (-El - a^) 2 ^ - Pix) 
dxi r\ 2 r\ 2 r\ 2 

by considering symmetry, we let 



(95) 



ri2 ■ (P2~Pl) _ P2-P1 _ P2-P\ _ 2(pi -p 2 ) 

V? i r 2 ~ r 2 o 2 — o 2 V yt) J 

'12 '12 J 'l2 J 'l2 



The second item in the left sideof Eq. (94) can be written as 



q\q2N 2 f , (Pl- P2) _ fi 2 • (p 2 - Pi) , 3 , 3 - 

F2 y 72 — ^ d°p 1 d°r 2 d°p 2 



AcmV 2 J m 1 12 



gig 2 ^v 2 /" , (pi -p 2 ) 2 ,3- ,3- .3 



6^F y / 2 ^f^^^rVV2 = -poa. (97) 
C7 s i can be called as the production of dissipative energy. Because 

Pl 2cm Pl rf 2 2cm r\ 2 P2 2cm r 2 2 y ' 

The third item in the left side of Eq.(88) can be written as 

/ h(F el ■ V gl +F e2 • V^U^frffxPfr = ' {Fe r2 ~ Fe2) - d 3 pid 3 f 2 d 3 p 2 

N 2 f - 

/ 2 V r - C/(r 12 ) • n 2 (F e i - F 32 )d 3 pid 3 f 2 d 3 p 2 = V • J s2 (99) 



4cm V^ 2 



It can be proved that Eq.(67) still hold when / is a vector, so let / = f(F e i — F e2 ), we have 

N2 f 1 J ^ f ^12^12 • < 2 [^ei(ri + (1 - A)f? 2) t) - F e2 (ri - \?{ 2 ,t)\ dU{r'{ 2 ) 



Js2 — 



[ iX S 



&cmV 2 J Q J r'{ 2 dr'{ 2 

x/ 2 (ri + (1 - X)f'{ 2 ,p 2 , fi - \^l 2 ,p 2 ,t)d 3 r'( 2 d 3 pid 3 p 2 

x/ 2 (n + (1 - AK 2 ,p 2 , ri - Afi 2 ,p 2 , t)d 3 ri' 2 d 3 pirf 3 P2 (100) 

By using Eq.(98) and (67) as well as relation Vf 1 J7(ri 2 ) = — V^ 2 U(ri 2 ), the fourth item in the left side of 
Eq.(87) is 

n 2 r 

-7^1 J ^i[V ?1 C/(r 12 ) • V fi + V r - 2 C/(r 12 ) • V^hd^d 3 ?^ 

= / ■ VrJJ(ri2)d 3 pid 3 r 2 d 3 p 2 = V • J s3 (101) 

cm v j r 12 



The transformation in the right side of Eq.(67) means that f\ — > fi + (1 — A)f^' 2 , r 2 — > fi — Ar^' 2 , so we have 
F12 = fi - r 2 — ► r*/ 2 and get 



fdX I fl2 „?f 2 ^"P Mn + (1 - Xy{ 2 ,P!,n - X^{ 2 ,p 2 ,t]d^p^p 2 

JO J r 12 r 12 " r 12 



83 2cmV 2 



= S^rI dX J ^^ i »+( i -^ft.^-^.A.*] rf, *? a d 'A rf, ft ( io2 ) 

By using Eq.(98) and relation Vf 1 f/(ri 2 ) = — V^ 2 ?7(r 12 ) again, The fifth item in the left side of Eq.(87) 
becomes 

N 3 f _ 

'W 2 J Ul ^^ U ^ ' V pi + Vr 2 U(r 23 ) • Vp 2 }.f3d 3 pid 3 r 2 d 3 p 2 d 3 r 3 d 3 p 3 

= f • Nr 1 U(r 13 ) - V F2 [/(r 23 )]dVid 3 r 2 d 3 p2d 3 r 3 dV 3 = V • (J s3 + J si ) (103) 

4cm i/g 7 r 12 

f gig 2 iV 2 f 1 y ^3^3 • (fi + (1 - A Mj - rj) dU(r'{ 3 ) 
8cmV 2 J J rg(fi + (l-A)f?3-f 2 ) 2 drf 3 

x/ 3 [n + (1 - A)fy 3 ,^i, r 2 ,p 2 ,fi - Af? 3 ,£», t]d 3 r'{ 3 d 3 p 1 d 3 p 2 d 3 r3d 3 p 3 (104) 

f _ gig 2 iV 2 f 1 f ?JA ■ [n - r 2 - (1 - A)f? 3 ] dU{r'j 3 ) 
85 8cmV 2 J 7 ^-^-(l-AJf^P dr'i, 

xf3[ri,Pi,r 2 + (1 - X)r 23 ,p 2 ,r 2 - Xr 2 ^,p 3 ,t]d 3 r 2 ^d 3 pid 3 p 2 d 3 r 3 d 3 p 3 (105) 

Similar to Eq.(99), the first item in the right side of Eq.(87) is 

N 2 r 

= -^Wl J h r{ 2 d Pld T2d P2 = V • Jse (1 ° 6) 

Similar to Eq.(98), we have 

?* = --£^vsf Q dX J ^VlV[Gi 2 (rl + (l-AK 2 ,t)-G 21 (rl-Art 2 ,i)]^^ 

xri' 2 / 2 (r! + (1 - Xy{ 2 ,p 2 , n - X7{ 2 ,p 2 , t^r'l^nd 3 ^ (107) 
The second item in the right side of Eq. (87) is 

n 2 r 

^2 j U X {K 12 ■ V ft + K 2 i ■ V p - 2 )/ 2 dVirfV 2 d 3 p 2 

- Sf / A *' - -v ■ /„ (108, 

x/ 2 (ri + (1 - AK 2 ,p 2 , fi - A^' 2 ,p 2 , i)d 3 r"d 3 pidV 2 (109) 

In the formula, we should let fi — > fi + (1 — A)f^ 2 , f 2 — > fi — Af^ 2 in the function Xi 2 and i? 2 i. Similarly, 
the third item is 

n 2 r - 

~WZ J Ul ^ Kl3 + F ^ ' + ( ^ 23 + F 2s) • ^p.]hd 3 pid 3 r 2 d 3 p 2 d 3 r 3 d 3 p 3 

qiq2N 2 f fi2 ■ r 12 {K 13 + F{ 3 - K 23 - F^ 3 ) 3 ^ 3 ^ 3 ^ 3 ^ 3 ^ f 

= 775- / J2 3 - j2 —d a pid a r 2 d' > p 2 d' > r 3 d a p 3 = V • J s8 (110) 

4cm l/g 7 r i2 



x/ 2 [ri + (1 - AK 2 ,p 2 , fi - Af5.2,p2j t]d 3 7{ 2 d 3 i?{d 3 p 2 (HI) 

We should also let fi — > ri + (1 — A)f^' 2 , r 2 — > ri — Ar^ 2 in ^13, i? 23 , ^13 and F 23 . For the last item in the 
right side of Eq.(87), we have 

n 3 r 

~2VJJ Ul[Vpl ' (G ' 13/3) + V?2 ' ( G '^h)]d 3 Pid 3 r 2 d 3 p 2 d 3 f 3 d 3 p 3 
= -f^£ / h ^ ' { f} 3 + ^ d 3 p ld 3 f 2 d 3 p 2 d 3 f 3 d 3 p 3 = V ft • J s9 (112) 

x/ 2 [ri + (1 - AK 2 ,p 2 , fi - \^ 2 J 2l t]d 3 ^ 2 d 3 p 1 d 3 p 2 (113) 
Let J s = J s i + J s2 + J S 3 + ^4 + Js5 + + + J.s& + J.s9, Eq. (87) can be written as at last 

d(p u s ) 



at 



V • ( Po u s V + J„) = p a s (114) 



Last let's discuss interaction caused by radiation damping forces. Suppose that the charges of particles 
arc distributed with spherical symmetry, we have 

U (r uPl ,t ) = — -2 (115) 

Similarly, producing eq.(50) by NU"(f\,pi,t)/2Vo and taking integral about <i 3 pi, we get 

1 1? / + |? / u 't - v ^ 3fi + w / t, " ,? " ' Vfi/,<i3fi 

TV 3 f a^ 2 7 

~2V? J U " S7 f 1 U(r 12 )-Vp 1 f 2 d 3 p 1 d 3 f 2 d 3 p 2 = -^ y ^"V fi • {G' 12 f 2 )dPp l dP?2<Pp2 

I U"(k 12 + F[ 2 ) ■ VpJ 2 d 3 p 1 d 3 f 2 d 3 p2 (116) 
Let energy density and energy fluid density of radiation damping force of unite mass are 

n 2 r 

U '=W*J U "*** (117) 



We get 



1 N 2 f -> 

P° u f V = 2 lp J VU"f 2 d 3 p ld 3 r 2 d 3 p 2 (119) 
V r - • (p u f V + J fl ) = V r - ■ {\{^) 2 J ^U"hd 3 n} 



lib f ( u "~ ■ v rJi + h- ■ v^t/Vpi (120) 

2 Vrf I m m 



Therefore, the second item on the left side of Eq.(116) can be written as 

n 2 r„„^ _ . ~ _ , -i -. , i jv 2 



2mVr 



/- 1 N z f 

U"pi • V^/idVi = V ft • (p u f V + J/i) - 2 0^72 _/ /iPi ' V r - 1 C/"rf 3 pi (121) 



Here 

2o 2 • 2o 2 ■ ■ 

V ft U"(n ,put) = V* (ai • pi) = ^-J- [pi x (V ?1 x a) + (pi • V Pl )ai] (122) 



Let 

g 2 A 2 
3c 2 m 2 F 2 



J fiPi ■ [Pi x (V?! x a) + (pi • V ?1 )ai]d 3 pi = PO071 (123) 
By using boundary condition, the second item in the right side of Eq.(116) can be written as 
A 2 f n 2 N 2 ( 

^72 J hFei ■ V^C/Vpi = 3 J m2y2 J /i^i ' Pi x ( V pi x S) + (Pi ' V ?1 )ai]d 3 pi - (124) 

The forth item on the left side of Eq.(116) can also be written as 

n 2 N 2 f 

V • J f2 = 3c 2 m2V r2 ) ./'2V r -C/(r 12 ) • V^W'dp^r^pz 

a 2 N 2 f ■ 
= 3 c2 m 2 V 2 j hVrMrn) • [pi x (V fi x S) + (pi • Vp-JaiJrfVi^f,^ (125) 

^ Z" 1 ^ /" ^3 ^12^12 ' [Pi x (V ft x 5) + (pi • V ft )5i] d^(ri' 2 ) 



J f^~6c^J dX J d ^ 2 ^ 
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x/ 2 [ri + (1 - A)ri' 2 ,pi,ri - \^ 2 ,p 2 ,t]d A ^ 2 d A p l d i p 2 (126) 
Similarly, the first and second items on the right side of Eq.(116) can be written as 

A 2 f 

P° a f3 = ^72 j hGi2 ■ V p - 1 ?7"d 3 pid 3 r 2 d 3 p 2 

72 J hGii ■ [Pi x (V ft x 3) + (pi • V # J3i]d 3 pld 3 f 2 d 3 p 2 = a /3 (127) 
N 3 f 

= y / 2 (#i2 + ^i 2 ) • V #1 i7"rf 3 p 1 d 3 f 2 rf 3 p 2 

a 2 A 3 /■ 

F2 j / 2 (^i 2 + F[ 2 ) • [pi x (V ft x 3) + (pi • Vpja^p^d 3 ^ (128) 



A : 



3c 2 m 2 t /2 



3c 2 mV 2 

Let J/ = J/i + J/2, 07 = 071 + 072 + 073 + 074, the energy equivalent equation of radiation damping forces 
is 

^ll+W-( Po u f V + J)=p a f (129) 
At last, the energy transition equation of hydromechanics can be written as 

+ V • ( Po uV + J) = -P:VV + Po a (130) 

In the formula, u = Uk + u v + u s + u/, J = ,1k + J v + Js + Jf, P = Pk + P v + P s + -P/, = 0fc + 0« + 0/- 
How to obtain the transport parameters from Eq.(130) remains to be researched later. 

7. The definition of non-equivalent entropy of general systems 

As well known, though the thermodynamics of equilibrium states had matured since the last stage of 
nineteenth century, the general theory of the non-equilibrium systems has not yet been established up to 
now. The key to establish the general thermodynamic theory for the non-equilibrium systems is to define the 
correct non-equilibrium entropy. But this is still an unsolved problem now. Though some special theories for 
special non-equilibrium processes have defined their non-equilibrium entropies, for example, the irreversible 
thermodynamics based on the hypothesis of local equilibrium ( 6 ), the extended thermodynamics ( 7 ), the 
rational thermodynamics , as well as the famous Boltzmann's non-equilibrium entropy, none of them are 
with general significance. 



In the equilibrium thermodynamic, the first law of thermodynamic or the law of energy conservation is 



dQ = dE + PdV = dE-F-df (131) 

In the formula, Q is total heat, E is total internal energy, P is pressure, V is volume and F is external force. 
For a single component system, the entropy function S is defined as 

TdS = dE + PdV = dE-F-dr = dQ (132) 

Because the concept of absolute temperature is based on equilibrium states, Eq.(132) is only suitable for 
equilibrium states. According to irreversible thermodynamics based on the hypothesis of local equilibrium, a 
system can be divided into many small cells, and each cell can be regarded as a small equilibrium system in 
which local equilibrium temperature T(r, t) can still be defined. Let S m (j, i), E m (r, 1), Qm(r, t) and E m (f, t) 
represent entropy density, internal energy density, heat density and external force density individually, the 
total entropy S, total internal energy E, total heat Q and total external force F of a system can be written 
as 

S(t) = J S m (r,t)d 3 r E{t) = j E m (f,t)d 3 f (133) 

Q(t) = J Q m (r,t)d 3 r F(t) = J F m (r,t)d 3 r (134) 

The first law of thermodynamics and the definition of entropy can be written as 

dQ m = dE m - F m ■ df (135) 

TdS m = dQ m = dE m - F m ■ df (136) 

As we known that if the interactions between micro-particles are known, the ensemble probability functions 
can be obtained. If the ensemble probability function is known, the heat density or internal energy density 
is also known. In this case, the local equilibrium entropy depends on the local equilibrium temperature. In 
the current theory of local equilibrium thermodynamics, the function form of local equilibrium temperature 
T(r, t) can't be determined by theory, but we can decide it through measurement point by point in principle 
by means of a small enough thermometer. Then we can determine the function form non-equilibrium entropy 
according to Eq.(136). 

Because the concepts of temperature and entropy are defined based on the concept of equilibrium, when 
the hypothesis of local equilibrium can't hold, there are argument about whether or not the non-equilibrium 
temperature and non-equilibrium entropy can exit at present, thought the concepts of external energy, heat, 
and force are meaningful. It should be seen that the concept of entropy was put forward only for the 
purports to describe irreversibility of macro-system's evolution, so that the second law of thermodynamics 
can be expressed in the clear form of mathematics. So it is necessary for us to define the non-equilibrium 
entropy form this angle. The problem is how to define it, instead of whether or not it exists. For the 
definition of non-equilibrium entropy, two conditions should be satisfied. The first is that when a system 
reaches equilibrium state, the non-equilibrium entropy should identify with the equilibrium entropy defined 
in the equilibrium thermodynamics. The second is that the definition should satisfy the principle of entropy 
increment in the non-equilibrium processes. Conversely, if a function can satisfy these two conditions, it can 
be regarded as non-equilibrium entropy. 

On the other hand, entropy is not a physical quantity that can be measured directly. T the concept of 
entropy depends on other quantities by means of them entropy is constructed. If the concepts to construct 
entropy are meaningful, entropy is also meaningful. Entropy is also a extensive quantity with additivity. So 
we can define non-equilibrium entropy density, just as we can define energy density, heat density and forces 
density in equilibrium states without considering equilibrium concept. 

As for non-equilibrium temperature, the situation is different. Temperature is a quantity that can be 
measured directly. The concept of temperature depends on equilibrium closely. If equilibrium does not exist, 
there is no the concept of temperature, for we can not use thermometer to measure temperature in a non- 
equilibrium system in which the hypothesis of local equilibrium loses its meaning. Meanwhile, temperature 
is strength quantity. We can't define temperature density. The concept of temperature always reflects the 
whole nature of big enough system. So we can say that the concept of non-equilibrium temperature is 
meaningless while the concept of local equilibrium does not hold. 

In this way, we do not use the concept of non-equilibrium temperature in the general non-equilibrium 
systems. In order to define non-equilibrium entropy, similar to the local equilibrium entropy and by the 



simplest method, we introduce an unknown function R(x, t) and define non-equilibrium entropy by following 
relation for a single component system 



RdS m = dQ m = dE m - F m ■ dr (137) 

Because R is not temperature, it can't be determined by the point to point's measurement in system. It is 
proved blow that we can decide the form of R function by the method of statistical mechanics, and prove 
the entropy increment principle of general non-equivalent systems by the connection of statistical mechanics 
and thermodynamic. When systems reach equilibrium states, we have R = T = constant. So the function 
S m defined in Eq.(137) satisfies two conditions mentioned above and can be used as the definition of non- 
equilibrium entropy of general systems. It can also be proved at last that the forms of non-equilibrium entropy 
and the function R are not unique, so this is also a reason why R can not be regarded as non-equilibrium 
temperature. Therefore, according to Eq.(137), we have 

flWr,*) - S m0 (r,to) = = mt) dt dt = W(r,t)dt (138) 

Here W(f,t) = Q m (f,t)/R(f,t), Q m (r,t) — dQ m /dt + VQ m -dr/dt. It means that when the system evolves 
from the state Co at moment to, it reaches the state C\ at moment t\. By using Eq.(133), the non-equilibrium 
entropy's increase is 

- S (t ) =J J dQ R m { ^ ) d3f = J [Smi(r,ti) - S m0 (f,t )]d 3 f (139) 

Now let's discuss how to decide the form of function R. Let fi(x\,pi,t) = f(x,p, t), satisfies Eq.(50). In the 
6-dimension phase space, the non-equilibrium statistical entropy density is written as Sn(ji • • •fjv,pi ■■ -Pn, t). 
In the 6-dimension phase space, the non-equilibrium statistical entropy density is written as S\(xi,pi,t) = 
S p (x,p,t). We define their relation as 

S p (r,p,t) = / _1 J f N (n ■ ■ ■ r N ,Pi ■ ■ ■PN,t)S N (r 1 ■ ■ ■ r N ,pi ■ ■ -p N ,t)dn 2 ■ ■ ■ fl N (140) 
The total entropy is 

= //jv(ri • • • r N ,pi ■ ■ ■p N ,t)S N {fi ■ ■ ■ f N ,Pi ■ ■ -pN,t)dQ.\ ■ ■ ■ Q N 

f(r,p, t)S p (f,p, t)d 3 rd 3 p (141) 

Comparing Eq.(144) with (133), the relation between non-equilibrium thermodynamic entropy density S m 
and statistical entropy density S p is 

S m (r, t) = j f(f,p, t)S p (r,p, t)d 3 p (142) 

For non-equivalent states, let dS = dSi + dS e . dSi is the entropy generation in the inside of the system, dS e 
is the entropy fluid from the outside to the inside of the system with relations ^6) 

' W ' .l 3 r ^l = -h m . d t (143) 



dt J dt 
The equilibrium equation of entropy is 

^ +V .j ro =<7 m (144) 

Similar to the local equilibrium theory, for a single component system, the entropy fluid and entropy gener- 
ation can also be written as 



Jm = S m V + i J a m = J ■ V-^ - ill : W (145) 



in the formula, J = Jk + J v + Js + Jf is the hot fluid discussed above, II is viscosity stress tensor. By using 
Eq.(50) and (141), we get 

1r = 1 + /^xw- /(-*! ■ v,,f + 5„ K + /f w 
-/</ 



[_ Vf • (S p ^/) + /-£ • V r -S p + S P X + /^]<^V (146) 
m m at 



In the formula 



K = -F el • V fi /i -iV y (i?i 2 + F^ 12 + F/ 2 ) • V?J s d 3 r 2 d 3 p s -N j V ft • (G 12 f 2 )d 3 f 2 d 3 p 2 (147) 



Eq.(146) can be written as 



/ [-V f • (S p -£-/) + • V r -S p + S P K + /^]d 3 p (148) 
J m m at 



(149) 



dt 

Comparing it withEq.(144), we have 

Jm = I f-Sd 3 p a m = j (/^ • V r ~S p + KS P + f^)d 3 p 

J m J m at 

By using Eq.(145), we get 

S m V+^J= [ f^S p d 3 p (150) 
it J m 

1 1^ f t j. p ^ a i j.dS p ^ 3 . 



J • v i? ~ i? n : w = y (/ m ■ vpSp + KSp + f ~df )dp (151) 

Put Eq.(142) and J = J k + J v + J s + Jf into Eq.(150), we can write 

J(f,t) = J J q (r,p,t)d 3 p (152) 
Then replacing the integral about d 3 p, we get the form of the non-equilibrium statistical entropy 

-w (153) 

G{7,P,t) = \ f-~ m Zl -U?,M (154) 
J (p — mvy 



Putting Eq.(153) into Eq.(151), we get 



In the formula 



— + A - VR = BR (155) 



K _JJGpfp_ n : W + ■ + KG + /ff )rf 3 p 

A ~ mffGdip B ~ J /G(Pp ( ' 

On the other hand, the viscosity stress tensor in Eq.(145) can be written as ft — PSij + £ij, here P is 
pressure and it is known or can be calculated in principle if / is known. Sij is the viscosity tensor and 

can also be calculated if / is known. So ft can be regarded as known quantity. Therefore, as long as the 
probability function / is known by means of statistical physics, we can achieve the form of R function from 
Eq.(155) in principle. Thus the forms of non-equilibrium statistical entropy in the 6-dimention phase space 
can be determined by Eq.(153) and the form of non-equilibrium thermodynamic entropy can be determined 
by Eq.(142). The form of non-equilibrium statistical entropy in the 6N -dimention phase space can also be 
determined by Eq.(140) in principle. 



The equilibrium states are discussed now. If we define the states with R =constant as the equilibrium 
states, according to Eq.(155), the equilibrium condition is B = . From Eq.(156), by considering the fact 
that we have V ^constant when a system reaches equilibrium states, the equilibrium condition becomes 

(\f^.V P G + KG+&d 3 p = (157) 
J m dt 

This means <j m — comparing with Eq.(149) when R ^constant. 

On the other hand, if R has nothing to do with space coordinate r but relative to time t, we have 
R = R(t) = T(t), T is absolute temperature. Suppose Q m also has nothing to do with space coordinate r 
but relative to time t, i.e, Q m — Q m (t), because 

dQ m _ dQ m df _ dQ m 

dt ~ dt + ^ m ' dt~ dt 

Put it into Eq.(139), we get 

*<« -*-/?/ -I his -/?£*-/* (i59 » 

This is just equilibrium entropy in equilibrium thermodynamics. Therefore, the temperature can be the 
functions of time in the processes of equilibrium states, as long as the systems are uniform. The processes 
are just the so-called quasi-stationary processes, but this point has not been noted clearly in the current 
thermodynamics. It is useful for us to calculate entropy function by using this nature. 

Now let we prove the principle of entropy increment in the adiabatic processes. That is to prove the 
following relation for the non-equilibrium isolated systems 

— = [ I^dV>o 

dt J R dt 

In principle, if the distribution function is known, we can prove it directly by the method of statistical 
mechanics, but this is impossible at present because of difficulty in mathematics. So we prove it by the 
method of connecting thermodynamics and statistical mechanics. 

Firstly, because Q m (t) and R(t) are the continuous functions of time, so S(t) is also a continuous function 
of time. Next, suppose an isolated system is in an equilibrium state at time to- At time to + dt a disturbing 
force is acted on the system so that the system becomes a non-equilibrium states. Then the disturbing force 
is removed immediately and the system would evolutes in the isolated state and reaches another equilibrium 
state at time t n . It is impossible to have another equilibrium state between these two equilibrium states. 
If there exists the third equilibrium state, it means that the system would become non-equilibrium from 
equilibrium without disturbing force form outside. In this case, the second law of thermodynamics would be 
violated. Thus, what we should prove is that the entropy would never decrease during the non-equilibrium 
process between two equilibrium states. By using the method of reduction to absurdity, we will prove below 
that AS should be a monotonously increasing function during the whole time to — > t n . 

According to the theory of equilibrium thermodynamics, because the system is in the equilibrium states 
at initial and final states, we have AS = S(t n ) — S(t ) > 0. If AS does not increase monotonously, when AS 
changes from dS > to dS < or from dS < to dS > 0, it would appear a state with dS = or dS/dt = 
at a certain moment U with to < ti<t n . According to equilibrium theory, it means that the system must be 
in the equilibrium states at moment t{. However, this is impossible as mention above. On the other hand, 
because F m = for an isolated system, we have Q m = E m and get according to Eq.(138) 

There are two ways to make Eq.(161) possible. The first is E m = or E m ^constant. However, this is 
impossible for / ^ constant, so E m ^ constant in the non-equilibrium processes. The second way is that 
E m jR is not a single valued function so that for any boundary condition f = fi and f = ?2, we always have 

J -j^d 3 r = W(r 2 ,t)-W(r u t)=0 (162) 



(158) 



(160) 



However, this condition can't be satisfied generally, for the probability distribution function is a single 
valued function in general. So Eq.(161)and (162) can not hold in the non-equilibrium processes in general. 
Therefore, dS/dt ^ in non-equilibrium processes and non-equilibrium entropy S must be a monotone 
function. On other hand, because AS = S(t n ) — S(t ) > 0, so during the process from time to to t n , S must 
be a monotonously increasing function. Thus, the principle of non-equilibrium entropy increment is proved 
for the non-equilibrium processes. 

At last we discuss the uniqueness problem of the definition of non-equilibrium entropy. Because probabil- 
ity distribution function / is unique, the form of function R is unique, so the form of non-equilibrium entropy 
shown in (137) is also unique. On the other hand, if there exists other form's non-equilibrium entropies S' m , 
we can always define them as 

R'dS' m = dQ m + dY m (163) 

If Y m is an unknown function, the definition is meaningless for there are three unknown functions and we 
can't decide them by the method mentioned above. If Y m is the function of known functions, for example 
energy density, heat density and force density and so on, we can also determinate the forms of S' m and 
R' by the same method. If dY m = and R' = T in equilibrium states, the form of S' m would be the 
same as equilibrium entropy in equilibrium state. If Y m and R' are also the single valued functions and 
Y m /R' ^ in non-equilibrium states, we can also prove the principle of non-equilibrium entropy increment. 
If all conditions are satisfied, we can also consider S' m as the non-equilibrium entropy density. On the other 
hand, from Eq.(137) and (163) we have 

RdS = R'd(S' m - U m /R') (164) 

S m (r, t) = J X d(S> m -Y m /R') dt + ^ = J« (SL ^ + *g_ )dt + Sm0 (165) 

In this way, both S m and S' m can be regarded as non-equilibrium entropy densities connecting by relation 
above. So the definition of non-equilibrium entropy is not unique. This is just the reason why we can define 
different non-equilibrium entropies in the current non-equilibrium thermodynamics. Because R ^ R' , so R 
can't be regarded as non-equilibrium temperature, for non-equilibrium temperature as a measurable physical 
quantity should be unique if it exists. Because non-equilibrium entropy is directly immeasurable, so it can 
be non-unique. In this way, non-equilibrium entropy is not a state function and its increment can also be 
non-unique, thought equilibrium entropy is a state function and its increment is unique. 

Now we have completed the reform of classical statistical mechanics. By considering retarded elcctromag- 
nctical interaction, we can also introduce asymmetry of time reversal into quantum theory. This problem 
will be discussed later. 
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